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Electronic structure of stripe ordered I^-^Sr^NiCU is investigated. The system with x = | 
is insulator, in LSDA+U calculations, and shows charge and spin stripe, consistent with the ex- 
perimental results. Highly correlated system of x = \ is studied by using exact diagonalization 
of multi-orbital many body Hamiltonian derived from LDA calculations and including on-site and 
inter-site Coulomb interactions. The fluctuation of the residual spin on Ni 3+ (hole) site couples 
with the charge fluctuation between Ni 3+ and Ni 2+ states and this correlation lowers the total en- 
ergy. The resultant ground state is insulator with charge and spin stripe of the energy gap 0.9 eV, 
consistent with observed one. The on-site Coulomb interaction stabilizes integral valency of each 
Ni ion (Ni 3+ and Ni 2+ ), but does not induce the charge order. Two quantities, inter-site Coulomb 
interaction and anisotropy of hopping integrals, play an important role to form the charge and spin 
stripe order in a system of x — |. 



PACS numbers: 71.27.+a,71.45.Lr,75.30.Fv,72.80.Ga 

I. INTRODUCTION 

The stripe order of charge and spin has been found in 
several layered perovskitesi^i 3 - and organic conductors^ 
Both of them are of pseudo two-dimensional (2D) elec- 
tron system and have strong Coulomb interactions com- 
pared to hopping integrals. A single band system can 
be a model of organic conductors, while, in layered per- 
ovskites, a number of relevant orbitals depends on fill- 
ing, crystal field, exchange splitting etc. The inter-site 
Coulomb interaction is essential to the charge stripe in or- 
ganic conductors and, then, extended Hubbard model is 
adopted to explain the mechanism of stripe order therefS 
In layered perovskites, however, an origin of the stripe 
order is still controversial particularly in the perovskites 
other than cuprate, while the order in cuprate is attract- 
ing much attention in conjunction with marked suppres- 
sion of T c with hole doping | and rich physics. 

Nickel compound La2_ a; Sr a; Ni04 (LSNO) is a typical 
system of static stripe order of charge and spini^i 9 -^ It 
is an insulator with the total spin S = in wide range 
of Sr doping x (0 < x < 0.9) i This stability of in- 
sulating phase is quite different from the variety of the 
phases in cuprate case; cuprate changes to metal, insula- 
tor and superconductor depending on the hole concen- 
tration. LSNO at x = | shows the highest spin or- 
der temperature T S q ~ 200K, because stripe order of 
charge and spin is commensurate with lattice periodicity. 
The periods of charge and spin stripe are not generally 
commensurate to the lattice. Incommensurability e in- 
creases with increasing x and satisfies e > x in the region 
^>x>0, e<xin the region x > |. Here, incommen- 
surability e is defined as displacement of peak positions 
of super structure from reciprocal lattice points. Increase 
of e saturates in the region of ^ > x > h with the value 



e ~ 0.44. 9 In the region x > |, there exists commensu- 
rate charge ordered phase without magnetic order, called 
a checkerboard type charge order, between T$o = 80K 
and T C o = 480K. 10 Another experimental fact in LSNO 
is the dependence of the ratio between two lattice con- 
stants along with c- and a-axis (c/a-ratio) on x. The 
observed c/a-ratio has a maximum value at x = \, be- 
cause holes are first doped into the x 2 — y 2 -orbital in a 
region x < \ and, then, additional holes are doped into 
the 3z 2 — 1-orbital when x > ^J^- 

We present two issues on LSNO in this paper by using 
both LSDA+U method and the exact diagonalization of 
many body Hamiltonian. The first issue is that the inter- 
site Coulomb interaction is essential to static charge order 
in doped layered perovskites. Not only Hartree energy 
but also correlation energy due to the inter-site Coulomb 
interaction (beyond Hartree-Fock approximation) is im- 
portant. Because electron configurations fluctuate be- 
tween hole (Ni 3+ ) and non-hole (Ni 2+ ) states by hopping, 
the correlation energy is maximized at x — h . In this sit- 
uation, we need to diagonalize the many-body Hamilto- 
nian to know the true ground state at x = | . The second 
issue is that the spin stripe order in LSNO is attributed 
to the structure of multi-orbitals. A spin moment on a 
hole site is strongly correlated with surrounding spin mo- 
ments on non-hole sites in the multi-orbital system with 
fractional occupation, which is a essential difference from 
the case of single-orbital systems. 

The stripe order in layered perovskites is, in some 
cases, attributed to the Jahn- Teller (JT) distortion ! 12 : 13 
In the present system, however, this is not the case, be- 
cause the JT distortion is not consistent with the local 
symmetry of the observed stripe order of charge and spin 
in x = | LSNO. In the E g type JT distortion coupling 
with the e g orbitals, two oxygen atoms at opposite po- 
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sitions, centering Ni site, displace in opposite directions 
with each other. Consequently, the periodicity should be 
doubled along with -Ni-O-- ■ --0-M- line. It contradicts 
to the observed tripled structure. Distortion with unger- 
ade mode can be consistent with the observed order but 
increases the total energy. The spin structure at x = \ 
is inconsistent with the JT distortion, too. 

There is another candidate for the origin of stripe or- 
der in layered perovskites. That is long-ranged Coulomb 
interaction^ Also it is well known that the inter-site 
Coulomb interaction stabilizes the charge order in quar- 
terly filled (single-orbital) extended Hubbard modeli&i^ 
Thus, the inter-site Coulomb interaction can stabilize the 
electronic structure of insulator with static stripe order. 
In the electronic structure of LSNO, degeneracy of 3d 
e fl -orbitals causes two important parameters: splitting 
between 3z 2 — 1-orbital and x 2 — y 2 -orbital and the hop- 
ping integrals. It is reported that different values of these 
variables bring the system different orderi 13 ' 16 We should 
get reliable values of these parameters from the first prin- 
ciples electronic structure calculations. 

The present paper is organized as follows. The charge 
and spin stripe of LSNO (x = |) is discussed in Sec. ILTl 
based on LSDA+U calculation. We show in Sec. IIIII 
the problem that charge ordered solution is not found in 
LSDA+U calculations. Then we explain that the LSNO 
of x = I is highly correlated system and the correlation 
energy in LSDA+U method is not enough to stabilize 
the charge ordered solution as a ground state. In later 
Sections, exact diagonalization of many body Hamilto- 
nian is employed to investigate the electronic structure 
of LSNO (a; = |). The Hamiltonian derived from LDA 
calculation is explained in Sec. IIV1 Section [V] is devoted 
to discussion of the charge stripe order and charge corre- 
lation functions in the system of x = g . The excitation 
spectra and the energy gap of the system with x = h 
is discussed in Sec. IVII The spin stripe order of LSNO 
(x = |) is discussed in Sec. IVIII Finally Sec. IVIIIl is the 
conclusion. 



II. CHARGE AND SPIN STRIPE OF LSNO 
(x = |) BY USING LSDA+U METHOD 

LSDA+U metho d 17 ! 18 in conjunction with the linear 
muffin-tin orbital method with the atomic sphere approx- 
imatio n 19 i 20 includes on-site Coulomb and exchange in- 
teraction U, J with rotational invariant form . 18 ' 21 This 
on-site Coulomb term in LSDA+U Hamiltonian is called 
a "Hubbard correction term" . LSDA+U also includes 
the inter-site Coulomb interaction by means of Hartree 
energy. Therefore, LSDA+U method can explain the 
physics of any charge ordered system, if Hartree energy 
is enough to describe them. 

For low doping systems (x — 0, ^) where the correla- 
tion energy induced by the inter-site Coulomb interaction 
is small, the results of LSDA+U calculations can explain 
well experimental results. Details of the calculations are 
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FIG. 1: (color online) Projected density of states for each spin 
(up and down indicated by arrows) of electron in Las Sr i MO4 

by LSDA+U method for (a) Ni 2+ and (b) Ni 3+ . Energy ze- 
roth is fixed at the top of the occupied bands. Inset shows 
the location of each site. Local magnetic moment of each site 
is -1.46, -0.94, 1.51, -1.51, 0.93, 1.46^ B , respectively from 
the left to the right. 



as follows. The values of U are chosen to be 7.5 eV for 
Ni 2+ ions22 and to be 7.0 eV for Ni 3+ ions^ 3 - The value 
of J are chosen to be 0.88 eV. These values are con- 
sistent with those of constrained LDA*il In the present 
system, NiOg octahedra are elongated-^ with the direc- 
tion of c-axis and tilted^ The elongation splits the two 
e g -orbitals by A = 0.97 eV in LDA calculation, where 
we denote the energy splitting between 3z 2 — 1-orbital 
and x 2 — y 2 -orbital as A. On the other hand, the tilt of 
Ni06 octahedra does not change the electronic structure 
in both LDA and LSDA+U calculations; the localized 
magnetic moments on the Ni ion are 1.56 /ib (with tilt) 
and 1.54 \lb (without tilt), and the energy gap 3.66 eV 
(with tilt) and 3.73 eV (without tilt). Therefore, the tilt 
of NiOg octahedra is neglected hereafter. We, then, fix 
the total volume of unit cell and set the c/a-ratio to be 
the observed value 3.26. Crystallographic coordinates of 
atoms are fixed at those of x — and, with changing x, 
positions of all atoms are scaled. The magnetic unit cell 
at x = is a \[2 x \/2 x 1 supercell, and a m\/2 x V2 x 1 
supercell at x = — . 

At x = 0, calculated values of the band gap (3.73 eV) 
and local magnetic moment (1.54 /xb) agree well with the 
observed values of 4 eV— and 1.68 /xsri respectively. 
Figure [T] shows the projected density of states of LSNO 
at x — i . The position of this x 2 — y 2 -orbital on the Ni 3+ 
site shifts to energy region just above the Fermi energy. 
The calculated energy gap is 0.10 eV. The resultant spin 
structure is such that antiferromagnetic domains are sep- 
arated by hole stripe, and the magnetic moment localized 
on each site is -1.46, -0.94, 1.51, -1.51, 0.93, 1.46 
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respectively, from the left to the right. The center of the 
hole stripe is located on the Ni 3+ ions and no neighbor- 
ing Ni 3+ ion exists. The spins on two Ni 2+ ions sharing 
the same neighboring Ni 3+ ions are anti-parallel to each 
other (inset of Fig. [lja)). Introduction of the multiple 
Slater determinant decreases the spin moments on Ni 3+ 
ions, because the spin configuration with the opposite 
spin direction on Ni 3+ ions gives the same energy with 
the present spin configuration. This order is consistent 
with the experimentally observed stripe and should be as- 
signed to the real ground state. However, the calculated 
lowest energy state is of no charge order and the energy 
is lower by 0.5eV/cell than that of the real ground state 
in Fig. [TJ The spin ordered alignment of this calculated 
lowest energy state is different from that of the charge 
ordered state shown in the inset of FigfIJa). The local 
magnetic moment of the calculated lowest energy state 
is, on each Ni site in Fig. [T^a), 1.48, 1.33, 1.48, -1.48, 
— 1.33, — 1.48/xb, from the left to the light. There are two 
problems responsible for the energy increase of the real 
ground state with the charge order. One is the absence 
of correlation energy arising from the charge fluctuation 
between Ni 2+ and Ni 3+ configuration. The other is the 
fact that LSDA+U method often underestimates the cor- 
relation energy of antiferromagnetic bonds. 



III. ELECTRONIC STRUCTURE OF LSNO 
(x = |, 1) BY USING LSDA+U METHOD 

LaSrNi04 (x — 1) is observed to be paramagnetic 
metal, the occupation of the Ni ion distributes homo- 
geneously and each ion is Ni 3+ . On the contrary, calcu- 
lated ground state is antiferromagnetic metal and the 
projected density of states is shown in Fig. [21 This 
contradiction between the observation and the calcula- 
tion is due to unrealistic stabilization energy of the mag- 
netically ordered state in LSDA+U calculation of large 
on-site Coulomb interaction U, J, against the paramag- 
netic state. If we assume that the ground state is rep- 
resented by multiple Slater determinants, paramagnetic 
state would be represented as a linear combination of ran- 
dom spin configurations. Then spin polarization within 
atomic sphere lowers the on-site Coulomb and exchange 
interaction energy and correlation energy stabilize the 
paramagnetic metal phase. However, in the LSDA+U 
method, the ground state wavcfunction is represented by 
a single Slater determinant. Consequently, the param- 
agnetic state is only possible in the state where all ions 
have no spin polarization. This increases the total energy 
much, coupling with large U, J. A dominant component 
of the real ground state at x = 1 should be a linear 
combination of Slater determinants which have a single 
electron on each site and spins are not ordered. 

The LSDA+U calculation fails to present the param- 
agnetic metal phase as the ground state at x = 1, but 
the correlation effects cause correct change of the split- 
ting of two e g bands. Figure shows that the 3z 2 — 1 and 
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FIG. 2: (color online) Projected density of states for each 
spin (up and down indicated by arrows) of electron in anti- 
ferromagnetic metallic LSNO(a; = 1) by LSDA+U method. 
Energy zeroth is fixed at the Fermi energy. Resultant elec- 
tronic structure is antiferromagnetic metal (see text). The 
hole states are 3z 2 — 1-orbital in contrast to a; = | case. Lo- 
calized magnetic moment on each Ni site equals to 1.02 /is 



x 2 — y 2 bands locate above and below the Fermi energy, 
respectively. Comparing Figs. Q] and [21 one can see the 
change of the hole characters from the case of x = | to 
that of x = i. In the x = -| system, the hole is doped in 
x 2 — j/ 2 -orbital, while it is doped in 3z 2 — 1-orbital in the 
x = 1 system. This change of the hole character agrees 
with the experimental results described in Sec. U The 
origins of this change are the strong on-site Coulomb in- 
teraction U between two orbitals and the difference of the 
dispersion of two orbitals/bands. In the LDA calculation 
of the system of x = 0, two e g bands cross the Fermi 
energy. The 3z 2 — 1-band is narrow and lying on the en- 
ergy region —0.62 eV^0.69 eV, and the x 2 — y 2 -band is 
wide and lying on the energy region —0.86 eV~2.29 eV, 
where the energy zeroth is fixed at the Fermi energy. 
When small amount of hole is introduced, the hole is 
doped into x 2 — y 2 -band, because x 2 — ?/ 2 -band is av- 
eragely located higher than 3z 2 — 1-band. At the hole 
concentration x = 1 , orbital polarization is maximized if 
only one band is occupied and the other is push up over 
the Fermi energy. Because the lowest energy levels of e g - 
band are mainly consist of x 2 — y 2 -orbitals, the x 2 — y 2 
orbital is preferable to be occupied. Thus, the hole char- 
acter changes from x 2 — y 2 -orbital to 3z 2 — 1-orbital at 
x = 1. In LSDA+U calculation, however the spin and 
orbital polarization could not be completely maximized 
and the system becomes metal. The potential correction 
in Table Q] shows this relation between the hole character 
and the strong on-site Coulomb interaction U. Poten- 
tial corrections are the derivatives of Hubbard correction 
term of LSDA+U Hamiltonian with respect to each ele- 
ment of local occupation matrix (c^ff^j on each atom, 

where a and (3 denote atomic orbitals and a spin. In the 
present system, potential corrections are diagonal ma- 
trices due to symmetry and only diagonal elements are 
listed in Table Q] In such case, each diagonal elements 
shows the energy level shift of respective orbital induced 
by on-site Coulomb interaction. Potential correction of 
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FIG. 3: (color online) Projected density of states for each 
spin (up and down indicated by arrows) of electron in anti- 
ferromagnetic metallic LSNO(x = |) by LSDA+U method. 
Energy zeroth is fixed at the Fermi energy. Resultant elec- 
tronic structure is antiferromagnetic metal (see text). Local- 
ized magnetic moment on each Ni site equals to ±1.3 /is 



3z 2 - 1 | (i = |) nearly equals to that of x 2 — y 2 | 
(x = 1) and that of x 2 — y 2 f (x — -|) nearly equals 
to that of 3z 2 - 1 t = !)• Therefore, the on-site 
Coulomb interaction applied to the narrow 3z 2 — 1- and 
wide x 2 — y 2 -band induces the difference of the hole char- 
acter between the low concentration (x = 3) case and the 
high concentration (x — 1) case. 

In high-doped x — \ case, the real system has a in- 
commensurate stripe type charge order, as is described in 
Sec. |H The ground state of LSDA+U calculation shows, 
however, no charge order. Consequently, the electronic 
structure of the calculated system is metallic as shown in 
Fig.H 

The charge order of the system of x — \ would be ex- 
pected to be a commensurate checker-board type order, 
because of following three reasons. The first is that the 
periodic boundary condition forces the charge order to be 
commensurate to the lattice, and once we assume the in- 
commensurability e to be a rational value \ , the resultant 
charge order is checker-board type. The second is that 
in the temperature range Xgo < T < Tco, the checker- 
board type charge order is observed in the real system 
with x = i. The third reason is that the positive value 
of inter-site Coulomb interaction stabilizes the checker- 
board type charge order. We must notice that the inter- 
site Coulomb interaction is included in LSDA+U Hamil- 
tonian by means of Hartree energy. 

We could not find a converged charge ordered solution, 
even as a metastable state, by using LSDA+U calcula- 
tion. Localized magnetic moment of the four Ni ions in 
the unit cell equal to 1.29, 1.32, — 1.29, — 1.32^s, respec- 
tively. Therefore, all the Ni ions are of the same type, 
allowing the difference of local magnetic moment by 3%. 
Figure [3] shows the partial density of state of each orbitals 



of Ni ions. 

The Hartree energy may be insufficient to stabilizes 
the charge order. In the real system, the electron lo- 
calized on the Ni 2+ site has extending tail on Ni 3+ site 
and the center of localized electron can be fluctuating 
among Ni 2+ and surrounding Ni 3+ sites, in order to lower 
the correlation energy. Because the ground state wave- 
function is represented by a single Slater determinant in 
LSDA+U calculation, the fluctuation of the center of lo- 
calized electron is neglected, then the correlation energy 
is underestimated. Moreover, in the Hartree energy, the 
contribution of the tail to the charge on Ni 3+ site and 
the contribution of the head to the charge on Ni 2+ site 
are treated as if they were separate electrons. Such treat- 
ment gives an additional self-interaction and causes the 
increase of the estimated energy. Because the number of 
nearest neighbor (n.n.) pairs of Ni 2+ and Ni 3+ ions are 
maximized in the checker-board type charge ordered sys- 
tem, above two reasons become serious for the increase 
of estimated energy. 

Now the prescription for resolving these problems are 
follows. We need to preserve the antisymmetry for ex- 
changing of any pair of single electron wavefunctions 
in order to remove self-interaction and need to employ 
as many configurations of Slater determinants as possi- 
ble in order to correctly estimate the correlation energy. 
Thus we employ the exact diagonalization of many body 
Hamiltonian in following Sees. 
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TABLE I: Diagonal elements of potential correction (in unit 
of eV). Potential correction is a derivative of Hubbard correc- 
tion term in LSDA+U Hamiltonian with respect to respective 
elements of occupation matrix. That gives orbital dependent 
potential corresponding to fluctuations of orbital occupation. 



IV. DOUBLE ORBITAL EXTENDED 
HUBBARD HAMILTONIAN FOR LSNO (x = ±) 
BASED ON LDA CALCULATION 

To understand the ground state at x = ^ , we adopt the 

extended Hubbard Hamiltonian H of e g electrons on two- 
dimensional (2D) simple square lattice derived from re- 
sults of LSDA+U calculations, and diagonalize this many 
body Hamiltonian exactly by using the Lanczos method 
with the inverse-iteration method; 
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where the braces (• • •) denotes the summation over n.n. pairs and the symbol t Slater-Koster type hopping pa- 
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rameters. The on-site energy e^s are determined by 
LDA calculation at x = j2£ t dda = -0.543 eV, 
tddS = 0.058 eV for the n.n. pairs and \t' dda + \t' ddS = 
—0.018 eV, t' dd7T = —0.023 eV for the second n.n. pairs 
and A = e 3z 2_ 1 — e x 2_ y 2 = 0.97 eV. The energy zero-th 
is set at the midst of e 3z 2_ 1 and e x 2_ y 2. The matrix el- 
ements of the intra-atomic Coulomb interactions Uup-yS 
are represented as functions of averaged Coulomb and ex- 
change interaction U, J, by using the same expressions as 
in LSDA+U method.— The values of U, J are chosen to 
be 7.5 eV, 0.88 eV, respectively, on all sites. The value of 
the inter-site Coulomb interaction V is chosen as 0.5 eV, 
except explicit indication of the value of V. Later we will 
explain the reason of this choice, V — 0.5 eV, in detail. 
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FIG. 4: (color online) Two spin configuration at x = = in 
a \/8 x \/8 cell: (a) Neel order and (b) one of charge and 
spin stripe order among doubly degenerate ones. Up and 
down arrows in (a) and (b) denote spins and circles in (b) 
denote holes. In Neel order, spin correlation of four nearest 
neighbors has negative sign. Those of two second nearest 
neighbors and a third nearest neighbor have positive sign. In 
stripe order, two second nearest neighbors have opposite signs 
with each other, (c) The hopping anisotropy of second nearest 
neighbors. 



Because two diagonal directions are inequivalent in the 
real LSNO [x — i), we introduce a parameter 5 show- 
ing the anisotropy of the second n.n. (diagonal) hopping 
parameter t'. We add to or subtract from t' 2 _, 2 „2 „,2 
as t'x 2 - y 2 x 2 ~ y 2 ^ depending on the direction, as de- 
picted in Fig. EJc). We investigate the 5 dependence 
in the range of < 6 < 0.02 so that the sign of the 
second n.n. hopping parameter does not change. This 
anisotropy 5 reserves translational symmetry. As a re- 
sult, finite value of 5 does not induce charge order, at 
least in the present range of S. For example, the differ- 
ence between the charge-charge correlation of the ground 
state at V — eV under S — 0.02 and S — is at most 3%. 
However, once charge order is induced by other quantities 
(in the present case, by inter-site Coulomb interaction 
V), finite value of 5 changes essentially the spin order of 
the system, as is shown in Sec. I VIII Finite value of S 
neither changes essentially excitation spectra and energy 
gap. The width of energy gap reduced by a few percent 
at V = 0.5 eV. Therefore, we fix the value of 6 at 6 = 
in Sees. IIV1 IVl and [VH We will change the value of 8 in 
Sec. I VIII where the change of spin order induced by S is 
discussed. 

We diagonalize the many-body Hamiltonian of a sys- 



FIG. 5: (color online) V dependence of total energies and 
charge-charge correlations at x = 5 and 5 = 0, where V is 
inter-site Coulomb interaction. The state labeled "S" is singly 
degenerate while the other state labeled "D" is doubly degen- 
erate. Energy-zeroth is set at the ground state energy. State 
S is a ground state on the region OeV < V < 0.41eV, State D 
is a ground state on the region 0.41eV < V < leV. Suffix 1, 
2, 3, 7 denotes following sites: origin, nearest neighbor, sec- 
ond nearest neighbor and third nearest neighbor, respectively. 
These site indexes are depicted in Fig. [6|a) . 
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FIG. 6: (color online) (a) Site index of the present system at 
x = \ in a v8 x y8 cell, (b) A type of hole distribution of 
the Slater determinants with the complete checker-board type 
charge order (SDCCO in short), where closed and open circle 
denote non-hole and hole Ni site respectively. Another type of 
distribution is obtained by exchanging the hole site and non- 
hole site, (c) The occupation distribution of SDCCO in (b). 
A number in the open circle shows the occupation number of 
respective site. 



tern of twelve electrons on the planar y8 x v8 supercell, 
where each site has two e g -orbitals. A periodic boundary 
condition is imposed, avoiding the bunching of electron 
at corners, and causes commensurate checker-board type 
order in the exact diagonalization result instead of in- 
commensurate charge stripe in real LSNO (x — |) as is 
mentioned in Sec. IIIII The total S 2 and S z of this sys- 
tem is invariant, due to spherical symmetry of the spin- 
space. We then restrict the Hilbert space so as to total 
S z = 0, which reduces the dimension of Hamiltonian to 
(i 6 C 6 ) 2 = 64,128,064. 



V. CHARGE ORDER OF LSNO (x = |) 

At first, we discuss the ground state of LSNO at x = 
h. With increasing V, the ground state changes from 
singly degenerate state (State S) in the region OeV < V < 
0.41eV to doubly degenerate (State D) 0.41eV < V < 
leV. Smoothness of the connection among respective 
series of States S and D with respect to V is discussed in 
Appendix [A] 

Figure [S] shows the V dependence of the total ener- 
gies and the correlation functions of charge fluctuation 
(ShiSfij), where Sfii = n,; — (hi) and 6 is fixed at 0. At 
V = 0.5 eV, charge correlations (ShiShj) of State D equal 
to -0.135, 0.093, 0.096 for n.n. = (1,2)), second 
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FIG. 7: (color online) V dependence of the weight of Slater 
determinants with complete checker-board type charge order 
(SDCCO in short), where V is inter-site Coulomb interaction. 
The state labeled "S" and "D" are the same states as in Fig. [5] 
5 = 0. 



n.n. = (1,3)) and third n.n. = (1,7)) pairs 

respectively, and this corresponds to checkerboard type 
charge order shown in Fig. [6jb). The checker-board type 
charge order in State D exists even at V — OeV, though 
the amplitude is very small. The charge correlations for 
n.n., second n.n. and third n.n. pairs at V = OeV equal 
to —0.074, 0.007, 0.019, respectively. These values satisfy 
the checker-board type charge order: The positive corre- 
lation between second n.n pair nearly equals to that be- 
tween third n.n pair, and the correlation between the first 
n.n. pair is negative. The charge correlations in State S 
are homogeneous in comparison with those in State D. 
The charge correlation in State S at V = eV of first, sec- 
ond and third n.n. pairs equal to —0.046, —0.013, —0.056, 
respectively. Therefore, State S at V — eV does not 
have any charge order. The absence of charge order in 
State S in eV < V < 0.4 eV is attributed to a slow 
increase of charge correlations with increasing V in this 
region of V. 

Now we discuss three features in Fig. [SJ two are about 
correlation functions of State S and one about the energy 
difference between States S and D; (i) Rapid increase of 
the charge correlation of State S in the region 0.4 eV < 
V < 0.6 eV, (ii) merging charge correlation of State S 
into that of State D, in the region V > 0.6 eV and (iii) 
the constant energy difference between States S and D in 
the region V > 0.6 eV. All these features are related to 
respective wavefunctions of States S and D. 

Let us define the ideal complete checker-board type 
charge order whose Slater determinant wavefunction 
(SDCCO) should satisfy the following two conditions, (a) 
and (b), then we specify the characteristics of real wave- 
functions of States S and D by using the overlap with SD- 
CCO's: (a) every site has occupation one or two, and (b) 
all neighboring sites of a singly occupied site are doubly 
occupied sites and vice versa. A type of their distribu- 
tion of occupation is shown in Fig. [6] (c). Assuming both 
transfer integral and on-site exchange parameter J equal 
to zero, SDCCO's would be the ground states, where 
all the spin configurations are degenerate, because the 
numbers of neighboring hole (Ni 3+ ) and non-hole (Ni 2+ ) 
pairs are maximized. In an actual calculation, trans- 
fer integrals are finite and the SDCCO's are not eigen- 
state. The value of V determines how much SDCCO's 
are hybridized into respective eigenstates. The more V 
increase, the more hybridization of SDCCO's are prefer- 
able energetically. 

In the region of eV < V < 0.4 eV, the ground state is 
homogeneously extending State S since the kinetic energy 



is a source of gain of the total energy. On the other hand, 
in the region of 0.4 eV < V, the ground state changes to 
the charge ordered State D since the source of energy gain 
is the correlation energy due to the inter-site Coulomb in- 
teraction V. The dependence of the total weight of SD- 
CCO's is shown in Fig. [7) Comparing Figs. [5] and [71 the 
above features (i), (ii) and (iii) appear as follows. Since 
the wavefunction of the ground state State S is homoge- 
neous, the overlap with SDCCO is small in the region of 
eV < V < 0.4 eV. We discuss more detail about the 
weight of SDCCO's of State S at V = eV in Appendix 
|U Then, in the range 0.4 eV < V < 0.6 eV, the weight of 
SDCCO's rapidly increases. Because both States S and 
D are eigenstates of H and orthogonal with each other, 
the coefficients of SDCCO's in State S are different from 
those in State D. In the range 0.4 eV < V < 1 eV, States 
S and D have the same charge order but the spin order 
is quite different, which is discussed in Sec. I VIII Then, 
the difference of the total energies between States S and 
D is attributed to the difference of spin configurations. 
As is discussed in Sec. Q] Tso < Tco m the real system 
with x = i and, therefore, the energy scale of spin or- 
der is smaller than that of charge order. In fact, we can 
see a large energy difference between States S and D in 
the region of eV < V < 0.4 eV, since the charge order 
different between States S and D here. Then , the charac- 
teristics of the charge order are the same in States S and 
D in the range 0.6eV < V < leV (Feature (ii)), as seen 
in the charge correlation functions and the same weights 
of SDCCO's in States S and D, and the constant energy 
difference between States S and D (Feature (iii)). The 
observed high temperature state (T > Tso) in LSNO of 
x = | might be a mixture of charge ordered eigenstates 
with significant weight of SDCCO's, including the dou- 
bly degenerate States D and the singly degenerate State 
S. 

As mentioned in Sec.Hl the real ground state of LSNO 
(x = ^) shows charge stripe order, incommensurate to 
the lattice. And as mentioned in Sec. IIII1 the periodic 
boundary condition forces the charge order to be com- 
mensurate charge order of checker-board type. There- 
fore, the ground state of the Hamiltonian should have 
checker-board type charge order, and we choose V — 
0.5 eV from now on. The value V = 0.5 eV is not very 
unrealistic because a recent report shows V = 0.34 eV in 
GW approximation of LaMn03. 27 



VI. EXCITATION SPECTRA AND ENERGY 
GAP OF LSNO (x = §) 

We show calculated single-particle spectral functions of 
the system with x = | in Fig. [51 in order to see whether 
the spectra have energy gap and the system is insula- 
tor as the observed case. Inter-site Coulomb interaction 
is fixed at V — 0.5 eV as discussed in Sec. [V] Three 
fe-points are chosen; k = (0, 0)i, (f , f )± and (f,0)±, 
where a is the nearest Ni-Ni distance, the positions of 
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FIG. 8: (color online) Single-particle spectral function of 
8 sites, N — 12 (x — ^) system with S — and V = 
0.5 eV. Three panels are corresponding to respective fe-points, 



(1,0)1, (i 



(0,0), from the top to the bottom. Each 



value of fc corresponds to different boundary condition. Two 
arrows near 10 eV shows the position of the highest ioniza- 
tion level 9.4eV and the lowest affinity level 10.3eV in case 
of k — (0, 0)i. Spectral functions are approximated by the 
polynomial of 800 degree in H , where H denotes the Hamilto- 
nian. A small imaginary part (0.01 eV) is added to the energy 
u! for smearing 5-function peaks. 



Ni ions are (n,m)a, and the two translation vectors are 
(— 2,2)a and (2, 2)a. Each fc-point corresponds to peri- 
odic or antiperiodic boundary condition along respective 
translation vector. In order to obtain continuous spec- 
tra, we introduce an imaginary energy r\ of 0.01 eV and 
the value of rj determines the resolution of the spectra. 

rrn, , r , , . c ^ total width of spectra 
ihc value of n must satisfies n > *- , 

where to denotes the the highest power of H, because m 
is the upper limit of the number of (5-functions in the 
spectra. In the Lanczos method, the truncation error 
easily breaks down the orthogonality of Lanczos vectors. 
Because calculating spectral function may need the large 
to, we employ the shifted-COCG method,— which is a 
variant of the conjugate gradient (CG) method. It must 
be noticed that the shifted-COCG method is numerically 
stable. Because shifted-COCG method requires an ref- 
erence energy in the energy region of a peak of spec- 
tral weight, we need to know rough profile of spectral 
weight. Therefore, the shifted-COCG method should be 
used only when we need very fine profile of spectral func- 
tion. More details will be explained elsewhere. 

In the present case, the "total width of spectra" is 
roughly estimated as E(Ni 2+ ) - E(Ni 3+ ) ~ U = 7.5eV 
for ionization levels, and affinity levels, separately. Then 
we choose m = 800. We calculate the spectral function 
of to — 160 by Lanczos method, and then, calculate that 
of to = 800 by the shifted-COCG method. Actually, the 
width of spectra is wider than U due to the mixing of 
higher/lower occupation configurations. The spiky struc- 
ture of spectra is an artifact due to the choice of smaller 
r\ = 0.01 eV. Our choice of the value of r\ is intending 
to show the gap structure between ionization levels and 
affinity levels, at every fc-point. We fix the value of hop- 
ping anisotropy 6 = 0. The hopping anisotropy S does 
not change the shape of the spectral function in the range 
of < S < 0.02. The width of the gap at k = (0, 0)± in 
Fig.[5]is smallest among three fc-points, where the highest 
ionization level (HIL) is located at 9.4 eV and the lowest 
affinity level (LAL) 10.3 eV. Therefore, we conclude that 
this system is insulator with an energy gap of 0.9 eV. 

It is important to know the symmetry of the single 
electron wavefunctions of HIL and LAL, in order to show 
that the "gap" is not an artifact of discretized fc-point. 



E b (eV) E t (eV) 



Peak 



The largest 
eigenvalue and 
its degeneracy 



Symmetry of 
eigenvectors 



HIL 
LAL 

HIL 
LAL 

HIL 
LAL 



at k = (0,0)i 



0.37 

2.87 
i 



9.15 9.55 
10.05 11.15 

at k = 
8.18 8.98 1.72 
10.02 11.02 2.99 

at k= (f,0)i 
8.18 9.18 1.88 
10.12 11.12 2.90 



0.09x4 
0.41x4 

0.63x2 
0.44x4 

0.38x4 
0.42x4 



\ 2 > 2 ) 

( — 1L\ 

\ 2 ' 2 I 



U' 4/ 

' —7r 3-7T x 

s 4 ' 4 * 



(l-o) 

(7,*) 



TABLE II: Bottom and top of the energy of peak-like struc- 
ture, area of the structure, the largest eigenvalue (and its de- 
generacy after "x" symbol) related to "single electron wave- 
function" (see text) and the representative k which shows the 
symmetry of the single electron wavefunction. In the right 
column, single k is written and the other equivalent (k) are 
omitted. HIL and LAL stands for highest ionization level 
(HIL) and lowest affinity level (LAL) respectively. Each fc- 
points corresponding to the boundary conditions are written 
in the leading lines. 8 = 0. 



The symmetry of respective single electron wavefunction 
can be labeled by crystal momentum fc. If the "gap" is an 
artifact and the ground state is metal, HIL and LAL are 
labeled by different k. Actually many levels overlap and 
form a single peak-like structure because of complicated 
interactions, especially in the energy range near LAL. 
Even the definition of the "single electron wavefunction" 
are not clear in such a circumstance. Therefore, we define 
a "single electron wavefunction" as follows, here. First 
we choose the bottom (2?b) and top (E t ) of a single peak- 
like structure. Next we integrate the Green function of 
matrix form in the energy range E\, < E < E t : 



1 f Et 

- / Im G^ j73 (w)du;. 



(2) 



and then we diagonalize it. This matrix corresponds to 



energetically partitioned \& ia Cjp) [Eb < E t < Ep) or 



CiaCjp^ (Ep < Eb < E t ). Therefore, the eigenvec- 
tor of above matrix corresponds to the single electron 
wavefunction related to the peak-like structure and the 
eigenvalue corresponds to the occupation number of the 
single electron wavefunction. We summarize in Table |TT] 
the Eb,Et, eigenvalues and the k representing the sym- 
metry of eigenvectors for the peak-like structure near HIL 
and LAL. At k = (0,0)^, HIL and LAL single electron 
wavefunctions share the same symmetry. This can not 
occur in the metallic system. Therefore, we can conclude 
that the split between HIL and LAL is not an artifact 
and the system is an insulator. 

Now we show that the inter-site Coulomb interac- 
tion V induces the energy gap between occupied- and 
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FIG. 9: (color online) Spectral functions of 8 sites, N = 12 
(x = \) system are drawn at 8 = 0. (a)Each panel corre- 
sponds to V = 0.5, 0.2 and eV from the top to the bottom. 
In order to show the overlap between the ionization levels and 
the affinity levels at V = eV, they are drawn separately in 
the upper and the lower halves of respective panels. Spec- 
tral functions are approximated by the polynomial of H 160 , 
where H denotes the Hamiltonian. Imaginary energy 0.01 eV 
is added to the energy iv for smearing 5-function peaks. Thin 
vertical line shows the highest ionization levels 9.4 eV, 7.9 eV 
and 6.9 eV (at V = 0.5, 0.2 and eV). (b) Each orbital com- 
ponent (3z 2 — 1- or x 2 — y 2 - orbital) of the spectral function is 
drawn separately, near the highest ionization level. Thin ver- 
tical line shows the highest ionization levels, same as Fig. (a). 
Both the highest ionization level and the lowest affinity level 
consist of the x 2 — y 2 orbital an 
k = (§, See also TableHU 
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FIG. 10: (color online) Spin correlation function (Si • Sj) with 
the anisotropy parameter 5 of the second nearest hopping 
integrals. The value of V equals to 0.5 eV. Inset shows site 
location. 



unoccupied- states, not on-site Coulomb interaction U. 
Figure [5] (a), (b) shows the "spectral function" with re- 
spect to State D with V = 0.5, 0.2 and eV. Note that, 
since State D is not the ground state at V = 0.2 eV and 

V = eV, they (for V = 0.2 eV and V = 0.5 eV cases) are 
not satisfying the definition of spectral function. These 
spectral functions are approximated by rational function 
of the degree of H 160 by using Lanczos method. The 
top panel of Fig. [9] (a) and the bottom panel of Fig. [8] 
are the same, except that the latter uses higher degree 
of H so ° and the shifted-COCG. They are sharing char- 
acteristic peaks of HIL and LAL. Therefore, 160 degree 
is enough to discuss the characteristics of spectra. At 

V = 0.5 eV (top panel), HIL and LAL are (x 2 — y 2 )- and 
(3z 2 — l)-orbitals, respectively, and this agrees well with 
experimental observation and the spectrum in Fig. [3] At 

V = 0.2 eV (middle panel), HIL lowers to 7.9 eV and the 
energy gap becomes smaller than that at V = 0.5 eV. 
At V = eV (bottom panel), HIL lowers to 6.9 eV and 
two excitation peaks overlap with each other in HIL and 
LAL, i.e the energy gap vanishes. Thus we can conclude 
that the system with x — \ becomes insulator due to V, 
not U, and the inter-site Coulomb interaction V causes 
the energy stabilization and the opening the gap of the 
ground state State D. The gap width in the bulk limit is 
discussed in Appendix [Cl 



VII. SPIN ORDER OF LSNO (x = §) 

We discuss the spin stripe order induced by the 
anisotropy of hopping integral between the n.n. sites. 
Introduction of the anisotropy does not contradict to the 
local symmetry of observed charge and spin stripe. This 
anisotropy with 8 may couple with orthorhombic lattice 
distortion of T 2g symmetry. Coupling to the distortion 
reminds us of JT mechanism. However, the T 2g distor- 
tion is not a JT mode for the e g electron, and therefore, 



this mechanism is not of JT mechanism. 

The spin correlation function (Sj • Sj) at 8 = equals 
to —0.37, 0.12, —0.07 respectively when V = 0.5 eV and 
shows a lack of long range order. Equivalence of two di- 
agonal directions is an origin of the two-fold degeneracy, 
as is shown in the result of k = (f , f )^ in Sec. IVl 

Because large amount of holes are doped in x 2 — y 2 - 
orbital, the electron in x 2 — y 2 -orbital on Ni 2+ can have 
a long tail extending to Ni 3+ site. The hopping part of 
the total energy per electron is —0.30 eV and its abso- 
lute value is as twice large as that in the parent material 
(x = 0, N = 16), —0.15 eV. Following the localized elec- 
tron picture, increase of hopping energy indicates larger 
hybridization of x 2 — j/ 2 -orbital between Ni 2+ and Ni 3+ . 
Then, spin on Ni 2+ site has strong correlation with that 
on Ni 3+ site. Thus the spin correlation in x = \ system 
is more complicated than that in x = system where the 
localized spin picture is applicable. 

The present parameter set in H is another origin of the 
complexity of spin structure. If J S> A, the majority spin 
is preferable in order to lower the exchange energy when 
the electron visited the Ni 3+ site. In contrast, if A ^> J, 
the minority spin is preferable in order to lower on-site 
energy. The present values of J and A satisfy neither of 
two conditions. In addition to that, transfer integral and 
the inter-site Coulomb interaction V has the same energy 
scale. This causes the competition between the extended 
state and the localized single electron wavefunction. All 
these competitions may cause the instability against the 
perturbation. 

In the case of 8 ^ 0, the ground state energy de- 
creases with increasing 8 and can be expressed as (36.75— 
104. 3<5 2 ) eV. The spin correlation function is depicted in 
Fig. [TDJ The value of the correlation function between 
nearest neighbors, (Si • S 2 ) = —0.37 ~ —0.38, does not 
change. The local spin moment is(Sj 2 ) = 1.31. If two 
electrons localized on one site in a triplet state, four elec- 
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trons are localized on the four n.n. sites, and the total 
spin S of this totally six electron system equals to 1, then 
(Si • S2) = —0.375, which is almost the same as the cal- 
culated value of (Si • 82)- In the range of 6 > 0.01 eV, 
correlation functions of the two types of second neighbor 
pairs have opposite signs and the structure of the cal- 
culated spin stripe order is consistent with the observed 
one. This critical value d = 0.01 eV is less than 2% of the 
absolute value of the n.n. hopping tdda, and thus small 
anisotropy changes the spin structure drastically. 

The spin order is quite different between systems of 
the single-orbital and of the multi-orbitals. The former 
on 2D square lattice is the Neel order on \[2 x \[2 cell^ 
and the latter is the charge and spin stripe. Once we 
assume only x 2 — y 2 -orbital at each site and put 8 elec- 
trons on 16 sites, we get the value of (Si • S2) equal to 
—0.03 which is very small compared to the value —0.38 in 
LSNO^S, This reduction of (Si ■ S 2 ) in the single-orbital 
system is attributed to the absence of localized spin on 
a hole site and consistent with the Neel order mentioned 
above. However, a spin remains on a hole site in LSNO 
(x = i). Since two neighboring sites of a hole site have 
spins of opposite direction with each other, the spin state 
of a hole site is a linear combination of up- and down-spin 
due to symmetry. Because no term in the Hamiltonian 
flips a spin on a hole site, two antiparallel-spin electrons 
on two hole sites must exchange their positions through 
hopping process, in order to create such a linear combi- 
nation. This hopping process lowers correlation energy 
and the ground state must be represented by the multiple 
Slater determinant. In all hopping processes to create the 
multiple Slater determinants, that of exchanging a hole 
and a non-hole site is most probable, due to Large U , 
and the number of these pairs of sites is maximized at 
x = h , This causes insufficient description of the ground 
state by single determinant in LSNO at x = h. 

Finally, we stress that the above stripe spin correlation 
induced by the anisotropy of the transfer integrals is only 
seen under the checker-board type charge order. When 
the anisotropy 6 = 0.02 eV is introduced to State S (see 
Sec. E} in the range eV < V < 0.4 eV, the resultant 
state shows the spin order of Neel order type, which is 
also seen in S = case in the range eV < V < 0.4 eV. 
The Neel order of half-filled 3z 2 — 1-orbitals appears un- 
der the homogeneous charge distribution and the small S 
does not affect the spin order. Under the charge ordered 
condition, the single electron wavefunction is rather lo- 
calized, if not site-localized, then the transfer integrals 
are reduced effectively. Then the small S changes the 
spin order drastically into the spin stripe. 



VIII. CONCLUSION 

We discussed the charge and spin order of LSNO (x — 
|) by using LSDA+U method and those of LSNO (x = 
5) by using exact diagonalization of the double orbital 
extended Hubbard model derived from LDA calculations. 



In the exact diagonalization, charge and spin order is 
discussed by using charge-charge or spin-spin correlation 
functions. Excitation spectra of LSNO (x = |) are also 
calculated to show that the system is insulator. 

In conclusion, LSNO with x — ^ and ^ are both insu- 
lator with charge and spin stripe order. In both systems, 
diagonal hole stripes are separately located on Ni 3+ site 
in order to reduce hole-hole interaction energy induced 
by inter-site Coulomb interaction V. We discussed the 
important role of multi-orbitals and mixing of multiple 
Slater determinants especially in high-doped x = i sys- 
tem, where the charge order is induced by the correlation 
energy of the inter-site Coulomb interaction V. Charge 
order and the inter-site Coulomb interaction V are di- 
rectly related to the energy gap in the excitation spectra 
of the system with x = \ . Spin stripe occurs only under 
the condition of the existence of the charge order, with 
a help of anisotropy S in diagonal hopping. Thus the 
spin stripe is determined by the electronic structure with 
smaller energy scale than that of the charge stripe. This 
is consistent with the observation of Tco > ?so- Though 
the spin stripe is related to anisotropy, the mechanism is 
not of Jahn- Teller type unlike the usual ordering in tran- 
sition metal oxides, because anisotropy in LSNO does not 
couple to the Jahn- Teller mode. 

The stability of LNSO in different hole concentration x 
depends sensitively upon several physical quantities, e.g. 
tiaj/3, A, U, J, V and S, and we believe that we have suc- 
cessfully shown the new scope of the combination with 
the first principles electronic structure calculations and 
the many-electron theory. We have developed a very use- 
ful novel tool for the extremely large matrix of extended 
Hubbard Hamiltonian, the shifted-COCG method, which 
will be explain more details elsewhere. 

Lastly we comment on layered cuprates. The band 
gap of La2Cu04 (LCO) is 2 eV, narrower than La2Ni04, 
and hole doping makes systems metallic— Due to these 
facts, the screened Coulomb interaction in LCO be- 
comes smaller than in LSNO and the energy gain by 
hole hopping is more important in doped cuprates. 
This may be one of reasons why the hole stripe in 
Lai.48Ndo.4Sro.i2Cu04 runs in a direction along the n.n. 
pair and the hole concentration is one per two Cu sites 
in the stripe.— 
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APPENDIX A: METHOD TO GET SMOOTH 
EIGENVECTORS WITH RESPECT TO A 
PARAMETER IN HAMILTONIAN 

We use Lanczos method and the inverse-iteration 
method with CG method alternatively in order to ob- 
tain the exact eigenvectors. These methods work well 
for getting ground states, but careful treatment must be 
taken for solving two smooth eigenvectors with respect 
to V in Sec. |V] Because both methods are the energy 
minimization process, iterative application of them, with 
no care, gives the eigenvector of lower eigen energy than 
that of the aiming eigen energy smoothly connected from 
the solved eigen energy at adjacent values of V, though 
one chooses the solved eigenvector at adjacent value of V 
as a starting vector. When the two levels are nearly de- 
generate, this problem becomes seriously important. In 
the present case, at least three energy levels cross with 
changing V from eV to 0.5 eV, because at V = 0.5eV, 
the lowest three levels are 36.7549, 36.7551, 36.7773 in 
the unit of eV, (the first one is a level of State S and all 
the three levels are doubly degenerate). State S locates 
at the 36.7939 eV at V — 0.5 eV, higher than these three 
levels, while State S is the ground state at V = eV. 

The first technique to avoid this problem is to use small 
number of the dimension of the submatrix tridiagonalized 
by Lanczos method, since Lanczos method has stronger 
tendency to reduce trial eigenvalue than that of inverse- 
iteration with CG method. We tried successfully the di- 
mension in the range ~ 80 depending on the case. The 
second technique is to control the trial eigenvalue so that 
the residual is minimized and the trial eigenvalue does 
not jump, in the whole processes. Here, the residual is 
defined as ||(.ff — (trial energy)) (trial vector)||. Because 
increase of the residual is a sign of the transition into 
lower level, once it occurs, we dispose a new (and possi- 
bly lower) trial eigenvalue. 

Lastly, the inner product among resultant eigenvectors 
of respective values of the specific parameter (V in the 
present case) must be checked. The large value of inner 
product assures the smoothness of connection to the re- 
sultant eigenvectors. In the present case, the inner prod- 
uct between State S at V = eV and that at V = 0.5 eV 
equals to 0.54 and the inner product between State D 
at V = eV and that at V = 0.5 eV equals to 0.73. 
These values are large enough to assure the smoothness 
of these eigenvectors, because no other eigenvectors at 
the end point V — 0.5 eV can have larger inner product 
than above respective states. 



APPENDIX B: WEIGHT OF THE SLATER 
DETERMINANTS WITH COMPLETE CHARGE 
ORDERED V = 

One might have a question why State S has a finite 
weight 0.045 of SDCCO's at V = 0, though the state 
has no charge order. We conclude that the SDCCO's 



are not overweighted by following discussion. Because 
the energy scale of the charge order is larger than that 
of spin order in the present system, we neglect the spin 
configuration in each Slater determinant, and count only 
the charge degrees of freedom. Strong U inhibits the 
configurations with site occupation more than two and 
zero and there exist 8C4 = 70 types of the charge con- 
figurations, where the site occupations are restricted to 
be one or two. A number of charge configurations with 
complete checker-board type order equals to two. As- 
suming that all configurations have the same weight due 
to the absence of charge order, resultant weight of SD- 
CCO's equals to ^ ~ 0.029, which is consistent with 
above 0.045. 



APPENDIX C: GAP WIDTH IN THE BULK 
LIMIT 

First we discuss the value of the energy gap 0.9 eV in 
two way, applying two extreme approximations, and show 
how inter-site Coulomb interaction opens the energy gap. 
One of two approximations is such that HIL and LAL are 
approximated by single particle excitation related to an 
extended hole/electron, and the other is such that they 
are approximated by a single particle excitation of a site 
localized hole/electron. In the actual picture, excitations 
are described by dressed quasi-particle/hole and not by 
a single particle/hole, as is discussed later in the present 
section. 

We begin with the case that HIL/LAL is approxi- 
mated by single particle excitation of an extended state. 
Due to large splitting A = 0.97 eV between 3z 2 — 
1- and x 2 — y 2 -orbital, hole is doped into x 2 — y 2 - 
orbital. Therefore, low energy excitation mostly con- 
sists of x 2 — y 2 -orbital and we neglect 2>z 2 — 1-orbita! 
in the present paragraph. The single electron wave 
functions of HIL and LAL belong to the same symme- 
try group labeled k = (±§,±§)±. as in Table [HI By 
unitary transformation, four plain waves labeled k = 
(±~,±~)i change the forms into four checker-board 
type wavefunctions cos(^x) cos(^y), sin(^x) sin(-|y), 
cos(^x) sin(^y), sin^a;) cos(^y). These four wavefunc- 
tions give the two types of charge distribution. One 
type of them is shown in FigJS] (b) and the other 
is obtained with exchanging hole and non-hole sites. 
Two wavefunctions cos(^x) cos(^y) and sin(^x) sin(Sy) 
do not interact with on-site Coulomb interaction 
U, because (cos 2 (fir) cos 2 (fy))(sin 2 (f a:) sin 2 (§y)) = 
16 sin 2 (7ra;) sin 2 (7ry) = on each lattice point. In- 
stead, the interaction between them with inter- 
site Coulomb interaction V is roughly estimated 
as — x — x AV x n = — , where n = 

n n n ' 

(the number of sites whose amplitude is finite) = 4 in 
the present case. This is occupied- and unoccupied- 
splitting. Assuming cos(^x) cos(Sy) is occupied, then 
sm(jx) sm(^y) is unoccupied and vice versa. The 
same situation occurs for another pair of wavefunctions 
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cos(^a;) sin(f-y) and sin(^x) cos(^y). Thus, the sys- 
tem becomes insulator with gap roughly estimated as 
— 0.5 eV, except V — eV where all the occupied- 
and unoccupied- states are degenerate (See FigOJb)). 
Though the correlation induced by inter-site Coulomb 
interaction causes the splitting between occupied- and 
unoccupied- states (the energy gap), this mechanism is 
different from the normal Mott insulator, induced by on- 
site Coulomb interaction V <22 In addition to that, there 
is another difference between the current mechanism and 
that of Mott insulator; the wavefunction of HIL and LAL 
have checker-board type charge order, not sharing the 
center of charge distribution in the current mechanism, 
unlike they shares the center of charge distribution in 
Mott insulator. However, there is a problem with this 
estimation of the gap. Following the present estimation, 
energy gap vanishes in the bulk limit, because the gap is 
inversely proportional to the system size. 

In the next estimation, HIL /LAL excitation is approx- 
imated by a single particle excitation of site localized 
hole/electron. Neglecting hopping integrals, the ground 
state of the present Hamiltonian consists of four Ni 2+ 
(spin triplet, u 1 v 1 ) and four doublet Ni 3+ (spin doublet, 
u 1 ) located on respective sublattice, where u and v de- 
note 2>z 2 — 1- and x 2 — y 2 - orbital respectively. All these 
components are assigned to the Slater determinants with 
complete checker-board type charge order (SDCCO's) in 
Sec. El Then the LAL is x 2 - y 2 -orbital of Ni 3 + and 
the HIL is x 2 — y 2 -orbital of Ni 2+ (both have a major- 
ity spin). Consequently, the energy gap is estimated 



as {[£(Ni 2+ )] +&V - [£(Ni 3+ ) + W]} - {£(Ni 2+ ) - 
J B(Ni 3+ )} = W = 0.5eV, where E(-) denote the ground 
state energy of Ni ion for each ionization state. This is 
an exact solution even in the bulk limit . 

The result of the present calculation is in the inter- 
mediate region of above two estimations; the calculated 
energy gap 0.9 eV is greater than 0.5 eV in extreme cases 
by the extended HIL/LAL approximation and less than 
2.0 eV by the site localized HIL/LAL approximation. 
Therefore, in the calculated excitations (with respect to 
State D at V — 0.5 eV in Sec. [Vj) are well described 
by the wavefunction in intermediate region between site 
localized and extended states. And the split between 
occupied- and unoccupied- states occurs due to long- 
ranged (inter-site) Coulomb interaction, though they are 
not sharing the center and does not have large overlap. 

Finally, we discuss on the bulk limit of the energy gap 
in the present calculation. Taking the limit such that 
transfer integrals go to infinity, the system becomes para- 
magnetic metal. Taking the limit such that transfer in- 
tegrals go to zero the system becomes antiferromagnetic 
insulator. Therefore, there exists a critical transfer in- 
tegral where the bulk limit of the system changes from 
insulator (metal) to metal (insulator). The high weight 
of the Slater determinant with complete checker-board 
type charge order (> 0.4 at V = 0.5 eV, as in FiglTJ 
strongly suggests that the present choice of transfer inte- 
grals makes the bulk limit insulator, and that is consis- 
tent with the real LSNO (x = 
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